Spectral statistics near the quantum percolation threshold 
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The statistical properties of spectra of a three-dimensional quantum bond percolation system is 
studied in the vicinity of the metal insulator transition. In order to avoid the influence of small 
clusters, only regions of the spectra in which the density of states is rather smooth are analyzed. 
Using finite size scaling hypothesis, the critical quantum probability for bond occupation is found to 
be = 0.33 ±.01 while the critical exponent for the divergence of the localization length is estimated 
aa V = 1.35 ± .10. This later figure is consistent with the one found within the universality class of 
the standard Anderson model. 
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The present work is concerned with level statistics in an 
Anderson type quantum percolation model. More specif- 
ically, we consider a single particle in a three-dimensional 
lattice with binary distribution of bonds and analyze (nu- 
merically) the distribution P{s) of adjacent level spacings 
s for bond occupation probabihties close to the critical 
one (which marks the metal insulator transition). 

Level statistics in quantum systems and its relation to 
random matrix theories constitutes an iipportant tool 
for understanding the underlying physicsH. In particu- 
lar, correlations between energy eigenvalues of a single 
quantum particle interacting with random impurities in 
the diffusive regime are consistent with the predictions 
of Gaussian matrix ensemblesoB. Recently, it became 
clear that in the vicinity of a metal insulator transi- 
tion (provided it exists in guch systems) there is a dis- 
tinct kind of level statisticsQ. In this novel statistics, the 
critical exponent for the divergence of the localization 
length appeals, in numerous expressions for the various 
correlationsDilj. Hence, it is difficult to perceive a ran- 
dom matrix theory which adequately describe this critical 
statistics, although some progress has been recorded in 
this directioiO. One of the clearest indications for the 
existence of a different statistics in the neighborhood of 
the metal insulator transition is displayed in the behav- 
ior of the nearest level spacing distribution P{s), which, 
for large level spacing s, falls off slower than GaussianH. 
This is found to be the case both for the-Aaderson metal 
insulator transition in three dimensionstS'tj as well as for 
the Hall transition in two dimensional^. 

One of the motivations for studying level statistics in 
a quantum percolation model is related to the question 
of whether it belongs to the same uniKscsality class of 
the Anderson model with site disordeJlj^Ej. The answer 
to this question is by no means clear, despite the fact 
that quantum percolation can be regapied as a special 
variant of the general Anderson modelEII. For example, 
in some quantum percolation models, the value of the 



critical exponent v for the divergence of the localization 
length, as can be deduced from the transmission of the 
systeiHyiis-ifound to be smaller than that of the Anderson 
modellla't^l. Our analysis suggests that for a tight binding 
model the critical exponent (as can be deduced from the 
level statistics) for site disorder and that for quantum 
(bond) percolation are nearly identical. 

Another motivation (upon which wc will not elaborate in 
this work) concerns with the fractal nature of the wave 
function near the critical point. In particular, if the crit- 
ical quantum probability for bond occupation (denoted 
hereafter as Pq) is only slightly higher than the classical 
one (denoted hereafter as Pc) then the critical wave func- 
tions live on a fractal object, and the geometrical fractal 
dimension becomes relevant. 

Let us start by introducing the quantum percolation 
model and then explain how the nearest level spacing 
distribution is computed. Our calculations are based on 
a tight binding Hamiltonian, 



h.c), 



(1) 



where (ij) denotes nearest neighbors. The hopping ma- 
trix elements tij are independent random variables which 
assume the values 1 or with probabilities p and q = 1—p 
respectively. The underlying lattice is a three dimen- 
sional cube of length L with periodic boundary condi- 
tions. The missing bond probability q plays the role of 
disorder strength. For each realization k of bond oc- 
cupation probability, p, the above Hamiltonian is diag- 
onalizcd exactly, yielding a sequence of eigenvalues E^, 
n = 1,2, ...L^. This sequence is calculated for A'^ different 
realizations, where TV = 3000, 1400,750,450,300 for the 
corresponding different sample sizes L = 7, 9, 11, 13, 15. 
This corresponds to 10^ eigenvalues for each sample size. 

The average density of states (DOS) for L = 13 as a 
function of p is presented in Fig. |^. The most notice- 
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FIG. 2. The level spacing distribution for L = 13. One can 
see the transition from a GOE distribution (indicated by the 
thick full line) towards a Poisson distribution (indicated by 
the thick dashed line) as p decreases. 



FIG. 1. The DOS for L = 13 as function of energy for 
different bond occupation values. The connection between 
various small clusters and peaks in the DOS are indicated in 
the figure. In the inset, an enlargement of the region around 
is = is presented. 



able feature is the appearance of a series of sharp peaks 
in the average DOS which increase as p decreases. This 
feature was already noted in Ref. ^ where the DOS for 
a quantum percolating system was calculated using the 
Sturm sequence method. The origin of these peaks is the 
formation of small disconnected clusters of sites in the 
sample. For example, a single site with no connecting 
bonds to neighboring sites always contributes an eigen- 
value e = 0. The probability for such a site is equal to 
(1— p)^, therefore one expects a contribution of L^{l—p)^ 
eigenvalues-jequal to zero to the spectrum. This is in 
agreementcll with the observed height of the central peak 
in Fig. |l| (the bin size is 0.072) and with its variation as 
function of p. Another prominent feature is the appear- 
ance of a gap iH_.the DOS which depends on p around 
the central peakEj which may be seen in the inset of Fig. 
|l|. A cluster of two sites connected by a bound has a 
probability of p(l — p)^° to appear and contributes eigen- 
values e = ±1 to the spectra. Similarly, clusters of three 
sites contribute e = 0, ±a/2 and clusters of four sites con- 
tribute e = 0, 0, ±a/2 if all the sites are connected among 
themselves and e = ±(3 ± V5)/2 if only three bonds are 
present. It is interesting to note that gaps seem to de- 
velop also around these peaks. 

Here we face the question of how to study a spectrum 
for which some of the levels form degenerate clusters. 
Indeed, one can apply the various statistical measures 
of level statistics only if the density of states is smooth. 



Looking at Fig. ^ one may concentrate on three such 
regions centered around (I) E — ±0.4, (II) E — ±0.8, 
(III) E = ±1.2 (the spectrum for an odd L with periodic 
boundary conditions is not symmetric). In each region 
a fixed number of levels are taken (15, 31, 57, 95, 145 for 
L = 7,9,11,13,15) and the spectrums unfolded by the 
usual procedure, i.e., Xi^i = Xi + Si and Si = n{Ei^i — 
iJi)/(£'j_i_L„/2j+i ~ -£'i-[n/2j)- In the data presented here 
n = 13 is used, but no significant difference is seen for n — 
9. Within these guidelines, the distribution of adjacent 
level spacings for each region, sample size L and bond 
probability p is then calculated. 

A plot of P(s) as function of the bond occupation proba- 
bility for L = 13 is displayed in Fig. ^ It can be clearly 
seen that the expected transition from a Wigner like be- 
havior for large p to a Poisson behavior for small p is 
manifested. One should also note that all curves seem to 
intersect at s ~ 2, which reminds us of the situation for 
the Anderson transition with on-site disordeiu. As has 
been shown in Ref. ^ a very convenient way to obtain 
the mobility edge as well as the critical exponent of the 
transition v is to study the parameter 7(p, L) defined as 



P{s)ds - e-"" 



(2) 



which characterizes the transition from Wigner to Pois- 
son. Denoting by S,{p) the localization length, this func- 
tion is expected to show a scaling behavior "f{p, L) — 
f{L/^{p)) which in the vicinity of the critical quantum 
bond probability pq is expected to behave asQ 



j{p,L)^j{pg,L) + C 
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where C is a constant. In Fig. || curves of 7(p, L) for 
different sample sizes L are plotted for levels in the first 
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FIG. 3. The scaling function 7 as function of p for different 
sample sizes for levels around E = ±0.4. A clear convergence 
of all curves at pq 0.33 can be seen, as well as the expected 
change in the size dependence of 7. 
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FIG. 4. A fit of the numerical data around E — ±0.4 to 
the scaling function represented by the curve. 



energy domain. It is noticed that the curves cross at a 
single point at which the order of heights with respect to 
L is reversed. This is an indication for the existence of 
finite size one parameter scaling behavior. Similar situa- 
tion prevails also in regions II and III. 

Based on finite size one parameter scaling analysis, the 
procedure for calculating the critical bond probability, as 
well as the critical exponent goes as follows. The quantity 
"f{p, L) is calculated for many pairs (pi,Li). It is then 
considered as a certain scaling function f{x) of the scaling 
variable x — L^^'^ {p — Pq). For x ^ 00 the system is well 
inside the diffusive regime and hence f{x) — > 0. On the 
other hand For x — > — 00 the system is well inside the 
insulating regime and hence f{x) 1. Practically, it is 
useful to shift the variable x to y{x) = where a 

and 6 are respectively the minimum and maximum values 
assumed by x. Evidently, y(x) ranges between —1 and 
1. Then one expands f{x) in a series of Tschebicheff 
polynomials T„[?/(a;)] [n — 0, 1, 2, ...if ). Minimization 
of the set of differences \f{xi) — ^{pi,Li)\ results in the 
unknowns Pq, v and the expansion coefficients (namely, 
the scaling function itself). In all cases, it is sufficient to 
cut off the number of polynomials &t K — Yl. 

The following results are obtained: for region \ pq = 

0. 335 ± .005 and 1/ = 1.32 ± .08, for region II pq = 0.33 ± 
.005 and v = 1. 35±. 10 and for region III = 0.325±.005 
and I' — 1.35 ± .12. As a measure of the quality of the 
fit the numerical data and the fitted scaling function are 
plotted in Fig. |. It can be seen that as one might expect 
v is the same for all the three regions, while there is a 
small shift in pq as E increases. The value of Pq and its 
dependence on E is in perfect agreement with previous 
numerical studies of quantum percolation systemsE2lE2l. 
On the other hand, v is not consistent with the different 
values of the critical exponent obtained for those systems, 

1. e.,v = 0.38 in Ref. |l| and = 0.75 in Ref. [l|, but 



is remarkably clpse, to its value for the on-site disorder 
Anderson modelc3CI ly — 1.5 ± .1. 

Another quantity which is sensitive to the critical expo- 
nent v is the behavior of the tail of P{fi) at the transition 
point. According to Kravtsov et. aB ln(— ln(P(s))) — 
(2 — 7) ln(s) + const., where 7 = 1— This is not 
an accurate algorithm to calculate v since it depends on 
the behavior of P{s) at the tail of the distribution, for 
which the statistics is rather poor. It is important to 
note that in Ref. ^ Shklovskii et. al predict 7=1 even 
in the critical region with no dependence on /i which is 
supported by sopae recent numerical work on the on-site 
Anderson modeO. Nevertheless for the quantum perco- 
lation model we obtain 7 = 0.68± .16, which corresponds 
to I' = 1.04^^34''. A better measure for 7 is the number 
variance S2(iV) which should behave as Y?{N) cx N"' , 
at least for moderate values o£JV in which an ad|di|tional 
linear term recently predictedcJ is not significaniliil. The 
logarithm of the number variance hi{Y?{N)) versus In(iV) 
is plotted in Fig. ^. Two different regions for which a 
linear behavior is observed can be seen (i) 1.5 < iV < 12 
and (ii) 20 < N < 45. In between, a jump in T?'{N) 
can be seen which might be associated with some small 
cluster peak. A linear fit in (i) gives 7 = 0.74 ± .02 cor- 
responding to V ^ 1.28t;Jg and in (ii) 7 = 0.76 ± .03 
resulting in v = 1.39!^'^g. All the above estimations of 

V fall within the range obtained from the finite size one 
parameter scaling. 

Thus, based on the analysis of various statistical prop- 
erties of the quantum percolation spectra, the critical 
exponent in the well behaved regions of the spectra is 

V — 1.35 ± .10. This, at least for properties connected 
to the energy levels, seems to put the quantum percola- 
tion system in the same universality class as the usual 
on-site disorder Anderson model. The previous studies 
calculated v via the transmission of the system at ener- 
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FIG. 5. The logarithm of the number variance as function 
of the logarithm of the average number of states for the re- 
gion around E = ±0.8. Linear fits were performed for (i) 
1.5 < iV < 12 and (ii) 20 < iV < 45. For 12 < iV < 20 a non 
monotonous behavior of the number variance is seen. This 
behavior is probably connected to a small peak in the DOS 
at ~ 0.83. 



gies very close to E = 0. As can be seen in the inset 
of Fig. |l| the DOS has a very strong p dependence in 
that region. Thus a large part of the dependence of the 
transmission on p is probably due to the change in the 
DOS and not because of some changes in the localization 
properties. Therefore, it will be very interesting to ex- 
amine the transmission in regions of E where the DOS 
has only a weak dependence on p. 
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